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Introduction 


Pnenoraena  of  acoustic  reflection  on  plane  boundaries 
between  fluid  and  elastic  nedia  have  been  studied  in  the 
literature  both  theoretically  and  experimentally.1'2  A 
striking  feature  is  the  appearance  of  the  so-called  Rayleigh 
wave,  a surface  wave  generated  in  the  elastic  medium  at  a 
critical  angle  of  incidence,  and  propagating  with  a speed 
usually  below  and  close  to  the  shear  wave  speed  in  the  solid. 
Most  of  its  energy  is  concentrated  in  the  solid  near  the 
boundary,  but  part  of  it  leaks  into  the  fluid,  causing  an 
attenuation  (complex  wave  number)  in  the  propagation  direction 
along  the  surface. 

An  additional  type  of  surface  wave  was  shown  to  exist  by 
Stoneley.3  This  wave  has  a propagation  speed  close  to  that  of 
sound  in  the  fluid,  is  unattenuatsd  in  the  direction  of  propa- 
gation, and  most  of  its  energy  is  concentrated  in  the  fluid 
along  the  ooundary. 

In  addition  to  these  genuine  surface  waves,  there  also 
exist  two  types  of  lateral  waves  on  a flat  boundary  between 
solid  and  fluid,  which  propagate  with  the  compressional  and 
shear  speeds,  respectively,  of  bulk  waves  in  the  solid.1 

Scattering  from  solid  elastic  cylinders  has  been  investi- 
gated by  means  of  the  l.’atson-Sommerfeld  transformation  applied 
to  the  normal-mode  series.1'4  In  this  case,  one  finds  "creeping 
waves"  with  a surface  wave  type  behavior,  divided  into  two 
classes:  tnose  with  speeds  close  to  the  elastic  bulk  speeds 
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(Rayleigh  and  Whispering  Gallery  waves) , and  those  with  speeds 
close  to  the  sound  speed  in  tne  fluid  (Stoneley  and  Franz 
waves) . Previously,  the  Rayleigh  and  Stoneley-type  waves  (i.e., 
tnose  tending  toward  the  Rayleigh  and  Stoneley  waves  in  t le 
flat  limit)  were  studied  by  Grace  and  Goodman5  and  by  Lapin6  by 
analytic  methods,  while  numerical  discussions  of  the  Rayleigh 
and  Whispering  Gallery  modes  (higher  order  modes  which  arise 
because  of  the  curvature  of  tne  surface)  were  given  by  Doolittle, 
et  al.4  The  latter  authors  also  treated  the  Franz  waves,7  i.e., 
higher  order  modes  in  the  fluid  whicn  arise  because  of  the 
curvature  of  the  boundary  (and  which  also  exist  on  an  impene- 
trable surface) . 

In  the  present  v/ork,  we  shall  establish  the  connection 
between  creeping  rave  and  fiat  surface  wave  theory  by  investi- 
gating the  limit  of  acoustic  scattering  from  an  elastic  cylinder 
whose  radius  tends  to  infinity.  In  Chapter  I,  we  calculate 
the  behavior  of  the  circumferential  wave  modes  for  large  cylinder 
radii.  Accordingly,  the  appropriate  Debye-  or  Airy-type  asymp- 
totic expansions  for  the  cylinder  functions  are  used  to  solve 
the  secular  determinant  for  the  complex  surface  wave  numbers, 
numerical  results  for  the  Rayleigh,  Stoneley,  Franz,  and  v;hispering 
Gallery  rave  numbers  are  obtained  as  a function  of  fluid  wave 
number  times  cylinder  radius  for  a solid  aluminum  cylinder  in 
water,  and  in  one  case,  also  in  vacuum.  In  Chapter  II,  using 
the  Watson-Sommerfeld  transformation,  we  examine  the  behavior 
of  the  solution  for  the  problem  of  a cylinder  scattering  radia- 
tion from  a line  source  in  the  fluid  as  the  cylinder  radius  tends 
to  infinity.  In  this  limit,  the  residue  sums  corresponding  to 
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the  different  classes  of  circumferential  waves  found  in 
Chapter  I are  shown  to  tend  individually  toward  the  different 
types  of  surface  waves  found  on  the  flat  surface.  In  this 
way,  the  transition  of  creeping  wave  to  surface  wave  theory, 
as  the  scattering  object  tends  toward  a flat  surface,  is 
established. 


Chapter  I.  Surface  Wave  Modes  on  Elastic  Cylinders 


The  complex  wave  numbers  of  circumferential  waves  on 
an  elastic  cylinder  in  a fluid  are  obtained  as  the  roots 
of  a 3 x 3 determinant  which  may  be  derived  in  various  ways. 
One  way  consists  in  assuming  interior  and  exterior  solutions 
in  a form  describing  circumferential  propagation8  [with  a 
time  factor  exp(-i uft)  suppressed]: 

f = (ks) 

(la) 


/tfV'f  ^(ferr> 


and 


(lb) 


$ = fjv H?(kr)  (id 

where  , A are  the  elastic  potentials  and  !§  is  the  velocity 

potential  in  the  fluid;  k is  the  acoustic  wave  number  in  the 

fluid , and  k_  and  k_  are  the  wave  numbers  of  bulk  longitudinal 
Li  X 

(compressional)  and  transverse  (shear)  waves  in  the  solid, 
respectively.  Matching  boundary  conditions  on  the  cylinder 
surface5  then  leads  to  the  secular  determinant.  Alternately, 
when  the  problem  of  sound  scattering  by  a cylinder  is  solved 
using  the  Watson-Somraerfeld  transformation,4  the  same  deter- 
minant appears  in  the  denominator  of  the  scattered  field, 
and  its  zeros  give  pole-type  contributions  which  represent 
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circumferential  waves.  The  scattering  problem  is  discussed 
in  Chapter  II. 

The  determinant  in  question,  as  a function  of  V , is4 
given  by 


/'V 


_,u 

<Xy 

*vT' 

CKyX 

rtv* 

(2a) 


where 


V = V [jTy  (Xt)  — Xy  ^Ty 

co*y  TyC*T^ 

«F-  -XrT^vC^+^VM. 


Here,  a.  is  the  cylinder  radius;  x=ka.=  <*><*/<•  where /C  is 
the  sound  velocity  in  the  fluid;  Xl^- “ 
where  are  the  bulk  elastic  velocities 


/CL\.  &+y)/ ^ ^ =Wf*  ( 3 } 

that  depend  on  the  Lame  constants  2^l  and  on  the  density 
fa  of  the  cylinder  material,  while  is  the  density  of 

the  ambient  fluid.  The  primes  on  the  cylinder  functions 
denote  derivatives  with  respect  to  their  argument. 

Our  subsequent  analysis  shows  the  existence  of  different 
classes  of  zeros,  corresponding  to  physically  different  types 
of  surface  waves  which  have  been  classified  in  the  Introduction. 
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Each  complex  root  Qf  the  equation  D(v)=0)  inserted 

in  Eqs.  (1) , yields  a circumferential  wave  with  wave  number 
v/<x~,  phase  velocity  C?-to)ojvr  and  linear  attenuation  Xl/a. . 
Excitation  and  re-radiation  of  these  surface  waves  take  place 
at  a critical  angled  given  by  5iV\Q'=  ^/(Z • 

We  shall  be  concerned  with  the  case  of  large  cylinder 
radii,  or  large  values  of  the  parameters  v'  X (h  id)  Xl  and  XT  ♦ 

For  this  purpose,  it  will  be  necessary  to  utilize  asymptotic 
expansions  of  the  cylinder  functions,  which  assume  different 
forms  in  different  regions  of  the  complex  v-plane,  mainly 
separated  by  anti-Stokes  lines.  In  our  case,  the  appropriate 
division  of  the  complex  v-plane  is  shown  in  Figure  1.  Only 
zeros  in  the  first  quadrant  need  to  be  considered,  those  in 
the  second  and  fourth  quadrant  leading  to  exponentially 
increasing  waves,  and  those  in  the  third  quadrant  differing 
from  those  in  the  first  only  in  their  sense  of  circumnavigation. 
Regions  1-4  are  separated  by  the  anti-Stokes  lines  of  H** OO  > 

, and  K? XA  , respectively,  on  which  also  the  zeros 
of  these  functions  and  of  their  derivatives  are  located.  The 
circles,  with  radii  determined  by 

| V- XJ. / = <r <>■'/») > | v-  x| - frCx'A)  ( 4 ) 

( Xi  - Xi_  or  Xy  ) , define  regions  within  which  Airy-type 
asymptotic  expansions  are  more  accurate.  Outside,  Debye- 
type  expansions  are  appropriate,  of  these,  transition-type 
forms  must  be  used  near  their  corresponding  anti-Stokes  lines. 
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In  addition,  anti-Stokes  lines  for  ^ originate  at  Xcyr 
and  run  along  the  real  axis  to  the  left,  with  the  zeros  of 
TvCXst)  being  located  along  them. 

In  the  following,  we  shall  consider  physically  different 
types  of  surface  waves,  corresponding  to  different  types  of 
zeros  of  in  an  individual  fashion. 

IA.  The  Rayleigh  Zero 

In  the  limiting  case  of  the  cylinder  radius  Q--*00,  corre- 
sponding to  the  case  of  a flat  elastic  half-space  bounded  by 
a fluid,  the  speed  of  the  Rayleigh  wave  is  known  for 
practical  cases  to  lie  somewhat  below  the  speed  yC-r  of  the 
elastic  shear  bulk  wave.1  The  corresponding  zero  of  VCy), 
in  the  case  of  large  but  finite  radius  of  curvature,  will 
then  lie  to  the  right  of  Xr  in  Figure  1,  with  an  imaginary 
part  that  puts  it  above  the  real  axis.  In  this  case,  the 
appropriate  asymptotic  expansions  of  all  the  cylinder  functions 
appearing  in  D(v>  are  of  Debye  type;  cf.  Appendix  A. 

Equation  (2a)  when  set  equal  to  zero,  can  in  general  be 
rewritten  in  the  form 

xkOO=aa(»,cja-{xl.Xr  -P(xl£i  *2v*[i-xu£(xu^/ 

{^t  +a  j0«t3  * Mi  -x^[i -xr^} 

(5) 

where 

•f  <&>=  ^J(X^/jy(xd> 

gC<>=  Tr"(k)/^(Xi) 

- •*£/#- 1 , 


(6b) 
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h(x)=  . 


(6c) 


If  the  appropriate  Debye  expansions  of  Appendix  A are 
inserted,  we  obtain 

■f&WJ + [ ~(Jm] 

h</)~  <-(i-r9'  ~^pp^  [~(r-o_  + Wl 


(7c) 


where  T-^/x  both  being  /v/&”00>  Using  these 

expansions  and  calculating  Eq.  (5)  to  lowest  order  in  X;.^ 
writing  yields  the  well-known  generalized 

Rayleigh  equation  for  the  flat  half-space  bounded  by  a 
fluid, 1 


which  has  as  one  solution  the  (complex)  Rayleigh  wave 
number  kfc=c%*.  to  be  considered  in  this  section. 

The  behavior  of  the  Rayleigh  wave  number  on  a 
cylindrically  curved  surface  of  large  radius  of  curvature 
is  obtained  by  retaining  terms  of  order  xf1  in  Eq.  (5)  when 
inserting  Eqs.  (7)  into  it.  With  an  iterative  procedure 
that  starts  from  the  flat-limit  Rayleigh  wave  number 


fefc'  “ •+  &(aTa) 


(9) 
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the  result  for  is  found  as 

r 


-% 


*>  Off  w 


Aftr  kT  X*L 


+ 2 


(10a) 


where 


**;=  fc-*?)*  **  - (A^i) 

(l-L^t)  . This  agrees  with  an  expression  obtained  pre- 
viously by  Lapin6 • The  derivation  of  Eq.  (10a)  involved 
expansions  which  are  valid  under  the  conditions 

Of  these,  the  case  is  most  stringent  because  of 

the  proximity  of  to  kT  in  physical  cases  of  interest. 

For  the  case  of  a solid  aluminum  cylinder  in  water 

(c=1493  m/sec,  cL  = 6420  m/sec,  Cp  = 3040  m/sec,  p^  = 

3 3 

1 g/cm  , p2  = 2.7  g/cm  ),  we  calculated  the  flat  Rayleigh 
limit  by  solving  Eq.  (8)  numerically,  with  the  result 

kR/k-  O.  SV3  + O.  O/SS"  c. 


M'A 


(10b) 


(10c) 


(10d) 


(ID 
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Subsequently,  we  calculated  Gr.  from  Eq.  (10a)  , using  this 
value  for  The  results  are  shown  in  Figure  2,  where  the 

trajectory  of  as  a function  of  the  parameter  k0L  is 

plotted  as  crosses  in  the  complex  plane.  The  conditions  of 
Eq.  (10c)  would  indicate  that  the  most  reliable  results  are 
those  for  k<k±  70.  Also  shown  in  the  figure  are  points  of  the 
trajectory  obtained  by  Ugincius1*'9  who  used  convergent  expan- 
sions of  the  cylinder  functions  for  ta-i25,  and  Debye  expansions 
for  ko->25  for  a numerical  evaluation  of  the  roots  when  Eq.  (2a) 
is  set  equal  to  zero. 

The  present  results,  for  large  values  of  ka,  appear  to 
be  a natural  continuation  of  Ugincius'  zero  trajectory,  while 
the  present  low-lta  values  might  be  less  reliable  because  of 
the  approximations  used  as  mentioned.  Note  that  our  results 
were  obtained  as  an  expansion  away  from  the  Rayleigh  limit, 
while  those  of  Ugincius  tend  towards  it  without  having  assumed 
it  as  a limit. 

It  might  be  noted  that  in  the  limit  of  an  elastic  cylinder 
in  a vacuum  (*V*°  ) , the  flat  Rayleigh  wave  number  f*n  as  well  as 
the  correction  term  €*  become  purely  real,  indicating  no  radia- 
tive losses  of  the  Rayleigh  wave  in  this  case. 

IB.  The  Stoneley  Zero 

In  the  limit  of  the  flat  elastic  half-space,  the  other 
solution  of  Eq.  (8)  of  physical  interest  is  the  Stoneley  zero,3 
a real  root  corresponding  to  a speed  somewhat  less  than  the 
speed  of  sound  in  the  fluid.1  Accordingly,  in  order  to  obtain 
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an  extension  of  the  Stoneley  zero  to  the  case  of  a finitely 
curved  cylinder , the  transition-Debye  asymptotic  forms  (as 
outlined  in  Appendix  A)  were  used  in  the  left-hand  side  of 
Eq.  (5)  because  of  the  expected  vicinity  of  the  zero  to  the 
anti-Stokes  line  between  Regions  3 and  4 of  Figure  1.  Debye 
forms , rather  than  Airy  forms  which  are  more  appropriate  in 
the  circle  of  Region  7,  are  used  since  we  again  want  to  obtain 
the  Stoneley  zero  in  the  curved  case  by  expanding  about  the 
known  flat-limit  Stoneley  zero.  Sufficiently  close  to  zero 
curvature/  the  radius  of  the  Airy  circle  in  the  v/kxx.  plane 
("reduced  plane”)  becomes  small  enough  for  the  Stoneley  zero 
to  lie  outside  of  it.  (Even  inside  the  circle,  Debye  forms 
are  not  incorrect;  they  just  become  less  accurate) . 

Ordinary  (non-transition)  Debye  expansions  were  used  for 
the  other  cylinder  functions  in  Eq.  (5) , so  that  the  same 
expressions  as  in  Eqs.  (7a,  b)  were  inserted  in  Eq.  (5) . 
However,  the  transition-Debye  expansion 


(12) 

leads  to  the  expression 

-h  e (i3) 

now  to  be  used  in  Eq.  (5) . Retention  of  the  lowest-order 
subdominant  part  in  the  transition-Debye  form  yields  the 
exponential  in  Eq.  (13) . 
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Using  an  approach  analogous  to  that  for  the  Rayleigh 
zero,  we  find  the  Stoneley  wave  number  as 

/*S'  - + +• 

where  k5  is  the  flat-limit  wave  number,  and  where 


(14) 


■ -^h^[^-^+x5jcST 


ST 

LXxu 


r„J\ 


>*x; 


fffelV*. 

+ kW  fc>  wl 

k**£J  bw  “fe"ibs:**34 


with 


-3llia4  ^ 


WL  kJ*T 

3i 


(15a) 


(15b) 


This  expression  differs  from  a previous  one 
given  by  Lapin6  in  some  details.  It  is  valid  under 
conditions  similar  to  Eqs.  (10c,  d)  with  kt  replaced  by 
ks  ; the  most  stringent  one  is  the  analog  of  Eq.  (lOd) . 

It  is  important  to  note  that  the  only  imaginary  correc- 
tion term  to  the  (real)  Stoneley  wave  number  for  the  flat 
case  is  the  exponential  in  Eq.  (15a) , which  resulted  from 
retaining  the  subdominant  term  in  the  transition-Debye 
expansion  of  Physically,  it  corresponds  to  the  fact 
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that  on  a cylinder  the  Stoneley  wave  radiates  off  tangentially 
into  the  surrounding  fluid  while  on  a flat  surface,  its  wave 
number  is  real  and  it  cannot  radiate  off  any  energy. 

For  the  solid  aluminum  cylinder  in  water,  we  solved  Eq.  (8) 
numerically  and  obtained 


which  agrees  with  earlier  results.10  Inserting  into 

Eq.  (15a)  yields  the  points  presented  in  Figure  3,  plotted 


as  a parameter.  The  imaginary  scale  of  the  figure  is  greatly 
expanded  because  of  the  smallness  of  the  imaginary  part  of  the 
zero.  Due  to  the  mentioned  conditions  of  validity  of  the 
approximations,  the  points  in  the  horizontal  portion  of  the 
figure  may  not  be  numerically  reliable.  It  should  be  noted 
that  Lapin's6  formula  would  give  the  imaginary  parts  of  the 
Stoneley  zeros  only  half  as  large  as  given  by  us,  but  we  believe 
this  to  be  due  to  an  error  in  Lapin's  printed  expression. 

IC.  The  Franz  Zeros 

This  type  of  zeros  arises  due  to  the  finite  curvature  of 
the  surface,  and  therefore  exists  even  in  the  case  of  impene- 
trable objects.  The  corresponding  surface  waves  get  on  and 
off  the  surface  tangentially,1  and  they  are  no  longer  present 
in  the  limit  of  a flat  surface.  Therefore,  one  cannot  expand 
the  positions  of  the  zeros  about  the  flat  limit,  as  had  been 
done  in  the  preceding  sections.  Instead,  we  shall  expand  about 

the  known  positions  of  the  zeros  for  either  a soft  or  a rigid 


(16) 


as  a 


trajectory  in  the  complex  hsjk.  plane  with  values  of  ka. , 
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cylinder  with  finite  curvature.  The  latter  are  given  by  the 
complex  zeros  of  tCfo  or  , respectively,7  which  are 

located  along  the  anti-Stokes  line  between  Regions  3 and  4 in 


Figure  1.  Explicitly,  they  are11  (n=  1,  2, 


) : 


Vr 


+ s-c*1) 

where  are  the  zeros  of  the  Airy  function,  and  ^ those 
of  its  derivative: 


(17a) 


(17b) 


(17c) 


For  the  case  of  the  elastic  cylinder,  we  use  the 
ordinary  Debye  expansions  for  the  cylinder  functions  in 

of  Eqs . (6)  , but  the  Airy-type  asymptotic  expan- 
sions (Appendix  A)  for  The  latter  are  valid  inside  the 

circle  of  Region  7 in  Figure  1 and  had  also  been  used  to 
obtain  Eqs.  (17) , but  the  results  link  up  smoothly  with  those 
for  the  zeros  outside  the  circle,  in  the  transition  region 
between  Regions  3 and  4 where  transition  Debye  expansions  are 
used.12  In  fact,  in  the  reduced  v/ka.  plane,  the  radius  of 
the  circle  shrinks  with  increasing  ka  at  the  same  rate  at  which 
the  zeros  tend  towards  its  center. 

One  finds  in  this  way: 1 3 

3*  pax  5T 


(18) 


15 


with  p^  and  »a  defined  in  Appendix  A.  When  this  is  substituted 
in  Eq.  (5),  we  obtain  the  equation 


where 


0 ='$  ['  + ('M"y 


3x 

+ &Cx"^) 


(19a) 


(19b) 


and 

s=  ri  [-  (/- j T^X-ci-  tci) 

(19o) 

The  method  to  be  used  for  solving  Eq.  (19a)  will 
depend  on  the  magnitude  of  P.  Using  Eqs.  (A8)  and  the 
relations 

y=x+  b'Cx’k') 

we  obtain  from  Eq.  (19b) , to  lowest  order  in  the  quantity  X~'3: 


r~  d'(**<rf'/K) 

(20) 

where  . For  typical  fluid-solid  interfaces 

one  has  £Kl  0.1  , so  that  P4 1 for  it)3  (for  an 

aluminum  cylinder  in  water,  X^,— 6000)  . In  this  case,  we 
follow  a method  of  Streifer  and  Kodis13-15  and,  letting 
we  expand  ^(’f)  and  about  , defined  in 

Eq.  (17c)  as  the  zeros  of  fltOl)  • Subsequently , the  quantity 
ij-^His  expanded  in  powers  of  W^.-LFexp&Q/i),  and  we  obtain 
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from  Eq.  (19a) : 

<^Cv>=  '7new/%  ie  +J-Vjxj^3 

Using  Eq.  (A8d) , an  iterative  solution  of  this  equation 
gives  the  result  for  the  Franz  wave  numbers  (labeled  by 

"=<A-V- 

vfa/o. 

where 


KL-Qt-k})^  0-sT) 

1--&-  2x.  AL 
f>  s X 

+ Kr/k? 

t — H XL_ Xr/lL~  i- %l£(3.-*L-<x!lb 

X’u  Xt 

This  result  was  arrived  at  after  one  iteration,  assuming 
that  *rc4«r4l,  and  that 


(21) 

(22a) 


(22b) 

(22c) 

(22d) 

(22e) 

(22f) 


These  conditions  are 
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met  for  the  material  parameters  of  typical  fluid-solid  cases. 
Accordingly,  the  above  asymptotic  expression  for  constitutes 

-#/j 

a series  whose  terms  decrease  as  powers  of  X 1 . 

Note  that  in  Eq.  (22b) , the  first  three  terms  agree 
exactly  with  those  of  the  rigid  zeros,  Eq.  (17a),  and  that  the 
material  properties  enter  only  through  higher-order  terms. 

As  in  the  previous  cases  of  the  Rayleigh  and  Stoneley 
zeros,  the  binomial  expansions  used  in  obtaining  Eqs . (22a,  b) 


under  the  conditions 

x > * 

(23a) 

(23b) 

x ^ z i 

(23c) 

They  impose  lower  limits  on  x,  depending  on  the  order 
of  the  zero.  In  practice,  these  are  found  not  to  be  very 
stringent,  the  most  stringent  one  being  Eq.  (23a)  with 
t=T.  In  fact,  the  lower  limits  of  validity  were  found  to 
be  as  low  as  ka^3  for  FI,  and  e.g.  ka^60  for  F5,  for  an 
aluminum  cylinder  in  water. 

Equation  (22b)  will  be  a useful  expansion  for  a large 
range  of  x,  as  long  as  the  numerical  value  of  z is  such  that 
subsequent  terms  decrease.  However,  if  x becomes  very  large, 
x > x^  (including  the  flat  limit  x-»°o) , we  will  have  r>  i from 
Eq.  (20).  In  this  case,  we  expand  fo'Cv)  and  Ai(n-)  in  Eq.  (19a) 
about  >r  , i.e.,  about  the  zeros  of  defined  in  Eq.  (17c). 
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Then,  the  quantity  may  be  expanded  in  powers  of 

yielding 


/ 


Using  Eq.  (A8d)  and  iterating  again,  we  obtain  for  the 
Franz  wave  numbers  Eq.  (22a)  with 


(25) 


which  was  arrived  at  after  one  iteration,  assuming  4 j_ 

and  t-  m.  With  this  assumption,  Eq . (25)  again  represents 

-//a 

a series  whose  terms  decrease  as  powers  of  X / . The  first 

three  terms,  in  this  case,  agree  exactly  with  those  of  the 
soft  zeros,  Eq.  (17b),  and  the  material  properties  enter  only 
in  higher  order  terms.  Conditions  of  validity  are  now  Eqs. 

(23)  with  ^ replaced  by  ^ . 

Numerical  values  of  Eq.  (22b)  are  shown  in  Figures  4a 
and  4b  for  the  two  lowest  Franz  zeros,  FI  and  F2 , for  an  alum- 
inum cylinder  in  water,  as  compared  to  the  zeros  for  soft  and 
rigid  cylinders  in  water.  The  zeros  are  plotted  as  trajectories 
in  the  complex  plane  with  values  of  ka  as  a parameter. 

Only  those  soft-cylinder  zeros  are  shown  which  do  not  interfere 
with  the  rest  of  the  figure.  It  is  seen  that  in  the  range  of 
x presented  here,  the  elastic-cylinder  zeros  lie  close  to  those 
for  the  rigid  cylinder  as  expected.  For  increasing  mode  number, 
corresponding  elastic  and  rigid  zeros  move  closer  together. 
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ID.  The  Whispering  Gallery  Zeros 

This  type  of  zeros  also  arises  due  to  the  finite  curvature 
of  the  surface,  but  is  associated  with  the  material  properties 
of  the  elastic  solid,  and  therefore  does  not  exist  in  the  case 
of  impenetrable  objects.  Since  they  are  no  longer  present  in 
the  limit  of  a flat  surface,  we  cannot  expand  the  position  of 
the  zeros  about  the  flat  limit.  Instead,  we  find  the  positions 
of  the  longitudinal  and  transverse  Whispering  Gallery  zeros  by 
expanding  about  the  known  positions  of  the  zeros  of 
and  , respectively.  The  latter  are  located  on  the  real 

axis  in  Figure  1,  to  the  left  of  and  XT , respectively,  and 
are  given  by 1 ^ O' = *>**)■• -V* 


where  are  defined  in  Eq.  (17c). 

IDl . The  Transverse  Whispering  Gallery  Zeros 
In  this  case,  we  calculate  the  trajectories  for  the  zeros 
which  tend  toward  Xt  for  large  cylinder  radii.  The  region  of 
interest  in  Figure  1 is  Region  6,  where  we  use  the  ordinary 
Debye  expansions  of  Eqs.  (7)  for  > and  K0<)  , but 


(26) 


the  Airy-type  asymptotic  expansion  (Appendix  A)  for  : 


with 

in  terms  of  in  Eq.  (5)  using  Eq.  (6b) . Substit 

these  expansions  into  Eq.  (5) , we  obtain  the  equation 


Substituting 
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where 


[(T-f-rT 

It 

J 4-  1 

(28b) 


p ^xT 

Using  Eqs . (A8)  and  the  relations 

y=  xT  + 

dsi)"-w'*  ***>, 

-i/j 

we  obtain  from  Eq.  (28b) , to  lowest  order  in  the  quantity  *t  • 
/}  = #*/*). 

' ' (29) 

Therefore  f^.>lfor  all  Xj>  1,  and  letting  we  expand 

Ai'Ol)  and  /Hfy)  about  ^ . Then  the  quantity  is  expanded 
in  powers  of  yielding 


<^0cTv>-  >+  IT' + ( %l^)rT~ 3 +•  Kxt^3)  . 


(30) 


Using  Eq.  (A8d) , an  iterative  solution  of  this  equation 
gives  the  following  result  for  the  transverse  Whispering 
Gallery  wave  numbers  ): 

tv 


-WTvv\~  V»oT >*/< 


Ola) 


where 


+ V*r  (^  - Vf  ^ er<*T-) 


(31b) 
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with  *r  defined  in  Eq.  (22c)  and 


(31c) 


Old) 


This  result  was  arrived  at  after  one  iteration,  assuming 
that  0(uCe{ f^l.  The  first  three  terms  agree  exactly  with 
those  of  the  zeros  of  > Eg-  (26)  , and  the  material 

properties  enter  only  in  higher  order  terms. 

The  binomial  expansions  used  in  obtaining  Eqs.  (31a,  b) 
are  valid  under  the  conditions 


They  impose  lower  limits  on  X,  depending  on  the  order 
of  the  zero,  the  most  stringent  condition  being  Eq.  (32c) . 
For  an  aluminum  cylinder  in  water,  the  lower  limits  of 
validity  were  found  to  be  kav87  for  WT,1  and  e.g.  ka^547 
for  WT,5.  Numerical  values  of  Eqs.  (31)  for  this  case  are 
presented  in  Section  IE. 


(32a) 


(32b) 


(32c) 


22 


ID2 . The  Longitudinal  Whispering  Gallery  Zeros 
In  this  case,  we  calculate  the  trajectories  for  the 
zeros  which  tend  toward  Xu  for  large  cylinder  radii.  The 
region  of  interest  in  Figure  1 is  Region  5,  where  we  use 
the  ordinary  Debye  expansion  of  Eq.  (7c)  for  h(x),  and  the 
Airy-type  asymptotic  expansion  (Appendix  A)  for  f(Xc)! 

p >«:<$>  (33) 

with  and  w\  defined  in  Appendix  A.  The  method  used 

for  finding  the  longitudinal  zeros  involves  an  expansion 
about  the  zeros  of  3^.(Xu)  which  lie  on  the  anti-Stokes 
line  for  (cf.  Appendix  A and  Figure  1) , and,  there- 

fore, would  necessitate  the  use  of  the  transition  Debye 
expansion  (Appendix  A)  for  4C*t)' 

£ C*tV  - v^s_lJT  + fr(Xr). 

(34) 

However,  the  use  of  Eq.  (34)  renders  the  solution  of 
Eq.  (5)  intractable  analytically  because  of  the  presence 
of  an  overlapping  double  infinity  of  zeros.  One  group  arises 
from  the  Airy  function  of  and  the  other  from  the  tangent 
function  of  Xj  . The  latter  group  corresponds  to  the  transverse 
Whispering  Gallery  zeros  which  lie  outside  the  circle  of 
Region  6,  and  they  link  up  smoothly  with  our  results,  Eqs.  (31), 
for  the  zeros  inside  the  circle.  We,  therefore,  approximate 
in  Ea.  (34)  by  the  appropriate  non-transition  Debye 
expansion  (Appendix  A) : 

i 0-Tt/^  frCx?)  . 


(35) 
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We  will  determine  the  validity  of  this  approximation  by  a 
self-consistent  check  after  we  have  calculated  the  desired 
zero  positions. 

Writing  in  terms  of  using  Eq.  (6b)  and  sub- 

stituting these  expansions  into  Eq.  (5) , we  obtain  the  equation 

A'Ofi/faCf)*  J"l  (36a) 

Whe'6  p = ^-4— Q-T3)'^ 

L - ^TXi-ryXi-TrV'  l 

^>7  ^ ja/.,-rA\  (36b) 


-ixi+tm  + Mis- 

using Eqs.  (A8)  and  the  relations 

dx 'xlS 

(hrfaxd)  - - &/*>**  *"  McSy 

we  obtain  from  Eq.  (36b) , to  lowest  order  in  the  quantity  >£ 

* (37) 

Therefore,  for  y^>  lj  [^>1  for  all  Xl>L,  and  letting 

n-  we  expand  flcfy)  and  A'c(r^  about  ^ . Then  the  quantity 
is  expanded  in  powers  of  , yielding 

i-M-  > +ir'+fr(*)n~3  ^ - cr<x*) . (38) 

Using  Eq.  (A8d) , an  iterative  solution  of  this  equation 
gives  the  following  result  for  the  longitudinal  Whispering 
Gallery  wave  numbers  (r\- *>*>...)  : 


(39a) 
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where 


K>l,  - “ 


,h“  + i C^TU>- 

40  W 


+■  i fi*u*T  T t 


u.L+v^ 


(39b) 

with  defined  in  Eq.  (22c)  and 


UL= 


A 


K=*u 


(39c) 


(39d) 


This  result  was  arrived  at  after  one  iteration,  assuming 
that  . The  first  three  terms  agree  exactly  with 

those  of  the  zeros  of  Eq.  (26),  and  the  material 

properties  enter  only  in  higher  order  terms. 

The  binomial  expansions  used  in  obtaining  Eqs.  (39a,  b) 
are  valid  under  the  conditions 


* > 

KN£|  > 

(40a) 

| \3A 

^7  1^1  , 

(40b) 

■ \3k 


(40c) 
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They  impose  lower  limits  on  X , depending  on  the  order  of  the 
zero,  the  most  stringent  condition  being  Eq.  (40c) . For  an 
aluminum  cylinder  in  water,  the  lower  limits  of  validity  were 
found  to  be  ka^23  for  WL,1  and  e.g.  ka^^l  for  WL,  5. 

In  order  to  determine  the  validity  of  our  approximation 
for  we  substitute  our  result  for  , Eg.  (39b))  into 

the  exact  expression,  Eq.  (34),  and  into  the  approximate 
expression,  Eq.  (35)  , and  compare  the  two  for  the  case  of  an 
aluminum  cylinder  in  water.  We  find  that  the  approximation 
imposes  an  upper  limit  on  x which  increases  monotonically 
with  mode  number  A.  This  result  is  illustrated  in  Figure  5, 
where  we  have  plotted,  as  a function  of  rt,  the  value  of  ka  at 
which  the  absolute  value  of  the  relative  error  in  the  approx- 
imation begins  to  exceed  25  percent.  Numerical  values  of 
Eqs . (39)  are  presented  in  Section  IE. 

IE . Discussion  of  Results 

In  this  section,  we  present  a graphical  comparison,  for 
the  specific  example  of  a solid  aluminum  cylinder  in  water 
(and  in  one  case,  also  in  vacuum) , between  our  analytic  results 
for  the  various  zeros  and  the  results  of  Dickey,16  who  used 
numerical  methods  to  solve  for  the  roots  of  D(>')=0.  Dickey 
did  not  use  Eqs.  (7) , but  calculated  the  asymptotic  expansions 
of  the  cylinder  functions  directly  using  the  Airy-type  expan- 
sions of  Appendix  A or  Watson's  formulation1 7 of  the  Debye 
expansions.  It  is  found  that  the  two  methods  complement  each 
other  insofar  as  the  analytic  results  can  be  more  easily  carried 
to  very  high  values  of  ka  (where  both  methods  become  more 
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accurate,  but  where  the  numerical  trajectories  of  one  type 
of  zero  often  become  hard  to  determine  and  to  identify  among 
the  variety  of  other  zeros) , while  the  numerical  results 
retain  their  accuracy  down  to  lower  values  of  ka  than  the 
analytic  ones,  due  to  the  various  approximations  made  in  the 
latter  method.  For  the  case  of  the  Franz  zeros,  both  methods 
are  accurate  down  to  very  low  values  of  ka  and  the  agreement 
between  the  corresponding  results  is  excellent.  For  the 
Rayleigh,  Stoneley,  and  transverse  Whispering  Gallery  zeros, 
an  apparently  smooth  transition  is  obtained  from  the  numerical 
results  below  ka'vlOO  (below  ka^200  for  the  Whispering  Gallery 
zeros)  to  the  analytic  results  for  the  higher  values  of  ka 
up  to  ka->oo  . Numerical  results  for  the  longitudinal  Whispering 
Gallery  zeros  are  not  yet  available. 

Figure  6 presents  the  complex  trajectories  of  ke'/fe.  for 
the  Rayleigh  zero  as  a function  of  the  parameter  ka.  The 
circles  represent  the  numerical  and  the  crosses  the  analytic 
results  (the  latter  being  also  shown  in  Figure  2) . The  agree- 
ment becomes  better  as  ka  increases.  The  circles  appear  to 
tend  towards  the  calculated  flat  Rayleigh  limit  for  ka-**o  (square)  , 
which  had  been  used  as  the  anchor  point  for  the  analytic  cal- 
culation, but  which  has  no  connection  with  the  numerical  calcu- 
lation. 

In  Figures  7,  dispersion  curves  for  the  Rayleigh  zeros  are 
shown  for  aluminum  cylinders  in  water  (Figure  7a)  or  in  air 
and  vacuum  (Figure  7b)  . The  Rayleigh  wave  phase  velocity  /C*' 
is  plotted  in  Figure  7a  relative  to  the  sound  speed  in  water; 
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i.e.,  essentially  the  quantity  -C*'/c ■=  k/kg/f.  is  plotted,  where 
kft/r=  ReO***}  . The  results  (circles:  numerical  results; 

curve:  analytic  result)  tend  towards  the  flat  limit,  , 

as  ka-v®.  In  Figure  7b,  we  plot  the  numerical  Rayleigh  zeros 
(circles)  for  the  case  of  an  aluminum  cylinder  in  air  £c=330 
m/sec,  ^|=  0.00129  j/cw\3  ) and  the  analytic  results  (curve)  as 
well  as  some  previous  results  of  Viktorov8'18  (crosses)  for 
the  aluminum  cylinder  in  a vacuum.  Here,  the  values  of  A. /t'/c* 
(i.e.,  normalized  to  the  flat  Rayleigh  speed)  are  plotted  vs. 

The  flat  limit  for  aluminum-vacuum  used  here  was  taken18 
as  0. 933£t-2836  m/sec. 

Figure  8 presents  the  numerical  results  (solid  circles) 
and  the  analytic  results  (crosses)  for  the  Stoneley  wave  phase 
velocity,  >Cs</c  , approaching  the  flat  Stoneley  limit 
0.9975  (calculated  by  our  numerical  solution  of  the  flat 
Rayleigh  equation)  as  ka-*4*.  For  values  up  to  ka  = 100,  the 
agreement  is  not  as  close  as  for  the  case  of  the  Rayleigh  zero, 
but  as  pointed  out  earlier,  the  analytic  method  should  become 
valid  for  the  Stoneley  zero  only  at  relatively  higher  values 
of  ka  as  compared  to  the  Rayleigh  pole. 

The  first  five  Franz  zeros,  together  with  the  Stoneley 
zero  again,  are  shown  in  Figure  9;  here.  Figure  9a  presents 
dispersion  curves  of  >C.p/c,  and  £s’/>C  plotted  vs.  ka  and  Figure 
9b  shows  the  normalized  attenuations  or  imaginary  parts  of  the 
wave  numbers,  and  Vs'i/ha  , plotted  vs.  ka.  The  agree- 

ment between  the  solid  curves  (analytic  results)  and  the 
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circles  (numerical  results)  is  excellent,  reflecting  the 
increased  overlap  in  the  range  of  validity  of  the  two  methods 
for  this  case. 

Results  for  the  Stoneley  zero  (solid  circles  and  crosses) 
have  been  entered  in  Figures  9a  and  9b  also.  While  its  dis- 
persion curve  is  very  similar  to  that  of  one  of  the  lower 
Franz  zeros,  it  may  nevertheless  be  clearly  distinguished  from 
the  latter  by  its  much  lower  attenuation  as  seen  in  Figure  9b. 

In  Figure  10a,  the  dispersion  curves  of  «C«K>r/c  for  the 
first  five  transverse  Whispering  Gallery  zeros  and  for 

the  Rayleigh  zero  are  plotted  vs.  ka.  The  solid  curves  are 
the  numerical  results,  while  the  long  dashes  correspond  to 
the  analytic  results  for  which  the  previously  mentioned  con- 
ditions of  validity,  Eqs.  (10c)  and  (32c),  are  satisfied;  the 
short  dashes  represent  the  analytic  results  for  which  these 
conditions  are  not  satisfied.  There  is  a smooth  transition 
from  the  numerical  to  the  analytic  results,  with  the  value 
of  ka  for  which  they  link  up  increasing  with  mode  number,  as 
expected.  In  Figure  10b,  the  normalized  attenuations 
are  plotted  vs.  ka.  The  numerical  results  (solid  curves)  are 
shown  for  the  first  four  zeros,  while  the  analytic  results 
(long  dashes)  are  shown  for  the  first  two.  Although  the  agree 
ment  between  the  two  methods  is  not  as  good  here  as  it  was  in 
the  dispersion  curves,  there  is,  nevertheless,  a smooth  transi 
tion  from  the  extrapolated  numerical  (short  dashes)  to  the 
analytic  results  at  high  ka’s. 
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The  analytic  results  for  the  first  seven  longitudinal 
Whispering  Gallery  zeros  are  shown  in  Figure  11.  The  dis- 


Figure  lib.  For  each  mode,  the  region  of  greatest  accuracy 

is  shown  as  a solid  curve  whose  lower  limit  is  determined 

from  Eq.  (40c)  and  whose  upper  limit  is  ka  from  Figure  5. 

max 

Thus,  our  analytic  and  Dickey's  numerical  results, 
arrived  at  independently,  are  in  good  agreement  with  each 
other,  increasingly  so  at  high  values  of  ka  where  the  quantities 
tend  toward  their  expected  flat  limits. 


persion  curves  of  -Cwu/c  vs*  ka  are  ?iven  in  Figure  11a,  and 
the  normalized  attenuations 


Chapter  II.  The  Scattering  of  a Cylindrical  Wave 
by  a Large,  Solid  Elastic  Cylinder 


If  a cylindrical  wave,  emanating  from  an  infinite  line 
source  of  unit  strength  at  o)  in  the  fluid,  is  incident 

upon  a solid  elastic  cylinder  of  radius  a.  (Figure  12) , the 
total  acoustic  pressure  at  point  in  the  fluid  is1'19 

(with  a time  factor  exp (-uot)  suppressed] 


i-  JL  *vf  H^X 


o 


(41a) 


where 

eo-t  \ n>o  ^ 

B„=T>nH^)+^n 

(41b) 


o?H*<a 

-.Tl 

0(*\ 

Uv— 

o 

<* 

and  1^,  C^1-  , and«^r*'  (1=^3  ) are  defined  in  Eqs.  (2). 


(41c) 
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Application  of  the  Watson-Sommerfeld  transformation1'19 
to  the  series  solution  then  leads  to 


P=P*-P= 


where 


p 


Bv 


SirtDv' 


cosy(>N^-^L  Hj(k^ 


cos  v(/fr-0lv. 


(42a) 


(42b) 


(42c) 


and  the  contours  C ' and  C0  are  shown  in  Figure  13  and 
result  from  opening  up  the  original  contour  C of  the  Watson 
transformation.  The  contribution  of  the  "background 
integral"  has  been  shown  to  be  small20  and  will  be 
neglected.  The  contour  C0  surrounds  the  zeros  of  ~DV  (first- 
order  poles  of  the  integrand)  discussed  in  Chapter  I. 
Splitting  J>x  into  integrals  over  contours  C ^ and  Ca  (Figure 
14)  and  applying  Imai's  transformation1'19 


cos  vi 


LvfT  « 

e.  cosv<*) 


-ie  T : 


sin'H'v 


(43) 
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to  the  integral  over  serves  to  split  off  the  geometrical 
part  of  the  solution  (which  no  longer  has  l/sm'lfv  in  the 
integrand),  thus  yielding  residue  sums  which  converge  on  both 
the  insonified  and  shadow  sides  of  the  cylinder.  The  geometri- 
cal part  pj  can  be  evaluated  using  the  saddle  point  method 
(corresponding  to  the  far-field  approximation  ) / where 

the  saddle  point  contour  C$  is  drawn  in  Figure  14.  Also  shown 
is  the  saddle  point  (to  the  right  of  ka)  which  yields  the 
incident  wave  and  the  saddle  point  yx  which  yields  the  geo- 


metrically reflected  wave 


1/19 


and  separates  the  two  types  of 


residue  sums  and  (arising  from  the  integrals  around  con- 
tours Cx  and  C4  ) . Higher-order  saddle  points  yield  waves 
which  are  transmitted  through  the  cylinder2  0 . V.'e  then  have 


Pi+?»-+'Pl 


(44a) 


with 


P*= 


* Sivx'ttV*  *^yrA  '‘HA 

-W-r  K u"Vt  \ u<Vk  \ 

i n Hy^)  Hyk(K)> 

vn>vs 


(44b) 


'-S  Du  * 


(44c) 


(44d) 


where 


DV=^/SV 


(44e) 
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In  the  limit  of  infinite  radius,  it  can  be  shown  that  yields 
the  corresponding  geometrical  portions  of  the  field  (i.e., 
incident,  geometrically  reflected,  and  transmitted  waves)  for 
the  flat  elastic  half-space  (cf.  Appendix  B and  Brill20'21). 

We  are  concerned  here  with  the  residue  sums  and  , which 
yield  circumferential  waves,  in  the  limit  of  zero  curvature. 

We  first  examine  the  limiting  behavior  of  the  saddle  point 
>5,  since  its  position  determines  which  residue  sum  is  used. 

The  equation  which  Vs  satisfies  is:19 

COS~*  4.  cos"1  V*.  — 2.  COS~‘ v3  + <£) 

fir  kro  fe.a.  * (45) 

In  taking  the  limit  of  Eq.  (45)  , the  following  changes 
of  variable  are  used22  (and  will  be  useful  later  on): 


v'  s ka.  sin  8” 
.> 


S = 


(46a) 


(46b) 


(46c) 

where  & is  the  angle  of  incidence  on  the  flat  elastic 
half-space.  Thus , CL , r,  and  r0  tend  toward  infinity  while 
the  source-surface  and  receiver-surface  distances  remain 
constant.  Keeping  terms  in  Eq.  (45)  through  , we  then 

find  that: 


i in  ©■<  = 


(47) 
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When  we  make  the  associations 


S-*X> 


(48) 


it  is  clear  that  Eq.  (47)  is  exactly  the  equation  satisfied 
by  the  saddle  point  &0  for  the  flat  elastic  half-space  (cf. 
Appendix  B and  Fig.  Bl) . Thus,  the  saddle  point  which 
yields  the  geometrically  reflected  wave  in  the  cylindrical 


metrical  reflection  in  the  flat  case;  the  two  are  related  by 
the  transformation  Eq.  (46a)  between  the  V-  and  D’-planes.  We 


which,  in  the  limit,  maps  into  the  flat  case  shown  in  Figure 
B2  and  discussed  in  Appendix  B.  Thus,  we  consider  a source- 
surface-receiver  geometry  for  which,  in  the  flat  limit,  all 
the  surface  waves  (except  the  Stoneley)  contribute  to  the 
field  at  the  observation  point.  From  Figure  14,  it  is  clear 
that  residue  sum  pA  includes  the  Rayleigh  and  Whispering 
Gallery  poles,  while  residue  sum  includes  the  Stoneley  and 
Franz  poles. 

1 1 A . The  Residue  Sum  pa. 

Using  the  expansion1 


case  tends  toward  the  saddle  point  ©i  which  yields  the  geo- 


will investigate  the  case  where  (cf.  Figure  14) 


(49) 


we  can  rewrite  T*  as 


(50) 
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V7e  then  separate  into  the  sums  pa  p over  the  Franz  poles 
and  over  the  Stoneley  pole.  Using  the  Franz  pole 

positions  Vpn  of  Eq.  (25) , the  appropriate  Debye  expansions 
(cf . Appendix  A)  for  the  Hankel  functions  of  (p  and  ^ , and 
Eq.  (46b)  along  with  the  relations  (cf.  Figs.  15) 


we  find  the  following  asymptotic  expression  for  the  Franz 
waves : 


(50a) 


> f.=  cos-'a-, 


(50b) 


rf  ea:m',vfe"  b 


2X 


(51a) 


where 


(51b) 


k'L  -fT1 


(51c) 


and  the  arc  lengths  S and  S/  are  shown  in  Figures  15. 
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contributes  an  algebraic  (non-exponential) 


The  term 

factor  (e.g.,  see  Nussenzveig ' s calculation  for  the  soft 
sphere23).  Thus,  the  yik'1'  Franz  wave  gets  on  the  cylinder 
tangentially,  creeps  clockwise  as  in  Figure  15a  (or  counter- 
clockwise as  in  Figure  15b)  around  the  cylinder  [traversing 
the  arc  length  S (or  s1  ) ] with  speed  and  attenuation 

k;.,  and  gets  off  tangentially  after  mr\  circumnavigations . As 
a tends  toward  infinity,  k;S  and  kiS1  also  become  infinite, 
so  that  in  the  limit  of  zero  curvature,  the  Franz  waves  are 
exponentially  damped  out  and  never  reach  the  observation  point 
Thus,  they  do  not  contribute  to  the  field  in  the  flat  case. 

In  an  analogous  manner,  we  find  the  following  asymptotic 
expression  for  the  Stoneley  wave  (cf . Figures  15) : 


oo 


m=-o 


Xe‘ 


i>rs, 

1ST 


(52) 

where  we  have  assumed  that  the  fluid  is  slightly  lossy  so 
that  the  flat  Stoneley  wave  number  ks  is  complex  ( 1^=  k s,.+»-  ksi) 
Thus,  the  Stoneley  wave  is  also  exponentially  damped  out  and 
never  reaches  the  observer.  As  the  cylinder  radius  tends 
toward  infinity,  the  Stoneley  wave  contributes  to  the  field 
only  in  the  case  of  glancing  incidence  ( Slg^ka.)  , in  which  case 
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the  method  of  steepest  descent  must  be  modified  to  take  into 
account  the  effect  of  a pole  near  the  saddle  point  ^ (cf. 
Appendix  B) . 

I IB.  The  Residue  Sum  pi 

Using  the  expansion  of  Eq.  (49)  , we  can  rewrite  as 

T w L J 

Hyn0*-v^  W y^O^  > 

Rc, (53a) 

where 


■PlC^= 


(53b) 


ofj* 

* 73 
CKy 

(53c) 


and  «y  * and  (X^1  (c-1,2,3)  are  defined  in  Eqs.  (2).  We  note 
that  Eq.  (5) , which  we  solved  for  the  pole  positions,  corre- 
sponds exactly  to  set  equal  to  zero.  We  separate  "p* 

into  the  sums  j>±  over  the  Rayleigh  pole  and  Pt>L>  and 
over  the  longitudinal  and  transverse  Whispering  Gallery  poles, 
respectively. 

We  calculate  first.  Using  the  pole  position 


vft'=kA'a 


> 


(54) 


where  is  given  in  Eqs.  (9)  and  (10)  , the  appropriate 

Devye  expansions  (cf.  Appendix  A)  for  the  Hankel  functions 
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of  \ , Y'  , and  r;  , and  E q.  (45c)  , we  find  that 

a jl[^a-yft'cos-,JLfc_'in 

i-* 


nZCk<>'* 


(55a) 


(55b) 


where  (cf.  Appendix  B) 

%=(^-^A=fec«Ot. 


(55c) 


Since  4^6  set  equal  to  zero  and  calculated  to  lowest  order 
in  X corresponds  to  the  generalized  Rayleigh  Eq.  (8) , it  can 
be  shown  that,  using  Eq.  (46a) , 

f;  fiWUv».  \<M  > 

(56) 

where  D+  is  defined  in  Eq.  (B3c) . Combining  Eqs.  (55)  and 
(56),  and  using  Eq.  (46b)  along  with  tile  fact  that  asymptot- 
ically (cf.  Figure  15) 

Sin&y sm\=  > 

cosS^cos©^  = 

we  find  the  following  asymptotic  expression  for  the 
Rayleigh  i.'ave  _ 

dl  , /o ~_l  L s/t) 

+ h-aS^j)  ^a't-*v^Va/ 

J wv=D  > (59) 


(57a) 


(57b) 
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where  the  arc  lengths  and  are  shown  in  Figure  16. 

Thus,  the  Rayleigh  wave  is  excited  at  the  critical  angle 
creeps  clockwise  ( S^)  or  counterclockwise  ( S*) 
around  the  cylinder  with  speed  and  attenuation  Xwvfeft. 

and  radiates  off  the  cylinder  at  the  same  angle  after  m 
circumnavigations . As  a tends  toward  infinity,  also 
becones  infinite,  and  therefore  the  imaginary  part  of 
causes  the  wave  which  creeps  around  the  shadow  side  of  the 
cylinder  to  be  exponentially  damped  out,  so  that  it  never 
reaches  the  observer.  The  same  argument  holds  true  for  the 
multiple  circumnavigations  (nv+O)  and  therefore  the  only  wave 
which  reaches  the  observation  point  in  the  limit  of  infinite 
radius  is  the  one  which  traverses  the  finite  arc  length 
and  corresponds  exactly  to  the  Rayleigh  wave  for  the  flat 
elastic  half-space  [cf.  Eqs.  (B8) ] : 


Pl>R' 


/Vl 


(59) 


By  using  the  Vronskian  relation 

<»>_  „ ]_ 


2 h 


H%<A 

and  the  fact  that  (cf.  Eq.  (53b)] 
we  can  further  simplify  our  expression  for  f>± 

7*1=  - i I ; ■£  e;^<f 

V\  w>=o 


(60a) 


(60b) 


Ua/D-3vA^' 
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We  calculate  the  limiting  behavior  (x-»*o)  of  using  the 

pole  positions  of  Eqs.  (39)  even  though  the  method  used 

to  calculate  them  imposed  an  upper  limit  on  X which  increased 
with  mode  number  n(cf.  Chapter  ID2.  and  Figure  5).  The  just- 
ification for  this  procedure  will  be  seen  later  in  the  calcu- 
lation. For  we  use  the  pole  positions  of  Eqs.  (31)  . 

Then,  using  Eq.  (46c)  and  the  appropriate  Debye  expansions 
(cf.  Appendix  A)  for  the  Hankel  functions  of  X,  ¥*,  and  r#  , and 
keeping  one  more  order  of  accuracy  of  yn  in  the  phase  terms 
than  in  the  algebraic  factors,  we  find  that 


HtLQ*)  Hvm  Oh*)  ~ 


L^T^- 


■V*“S  ST-flJexp^R.R.-)^ 


(62a) 


Vn~  y (62b) 


where  (cf.  Appendix  B) 


X,T=  Cfe-feV)A=  k “S  \t  ■ 

In  order  to  calculate  the  quantity 


(62c) 


(63) 


we  use  the  following  asymptotic  expression  for  a* 
[i.e.,  right  hand  side  of  Eq.  (5)]: 


Ij<y> 


/■nJ 


j^T  [>+ 


+ yit 


(64) 
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where  — and  l=LjT  > is  defined  in  Eq.  (6a)  . 

The  following  Airy-type  asymptotic  expressions  for  $(xC) 
and  its  derivative  will  also  be  helpful  [cf.  Eqs.  (27)  and 
(33) ] : 

(65a) 


Plity  ® 


^ Aitya*  M * 


/* 


rsj  — 


'nj  — 


ty-l  Ai<£>J 


(65b) 


where  we  have  calculated  the  leading  order  behavior  of 
these  quantities.  Then  [ •f1Xyv\)  ]”'  for  is  found  by 

using  the  asymptotic  expansions  of  Chapter  ID2.  and  Eqs. 
(64)  and  (65)  along  with  the  fact  that  [ef.  Eqs.  (36)] 


■f  frc) 


_ l 

tw- 

(66) 


The  result  is 


V 

<*t 


fz.()-0  fa«L-*f)  K (67) 

Similarly,  [ ]”'  for  y„=vwTln  is  found  by  using  the 

asymptotic  expansions  of  Chapter  ID1.  and  Eqs.  (64)  and 
(65)  along  with  [cf.  Eqs.  (28)] 


and  the  result  is 

1-1 


JL 


(68) 


(69) 
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Finally,  using  Eqs.  (62)  and  (46b)  and  the  fact  that 

asymptotically  (cf.  Figure  17) 

8n  — • ft  — S©  — S 

> (70a) 

cos©/t/^  C05&r=  !*£-  = .§-> 

do  d (70b) 

(where  C^-  is  the  excitation  angle  of  the  A Whispering 

Gallery  wave)  we  find  the  following  asymptotic  expression 

for  the  longitudinal  and  transverse  Whispering  Gallery  waves : 

--  ^ 

Mvr  . 


fi  ,T~  - i-[^(v3"'{e; - Kt  sst]  sL; 
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W\-b 


where  ( v^T  ) ]'*  is  defined  in  Eqs.  (67)  and  (69), 
is  the  number  of  longitudinal  (transverse)  Whispering  Gallery 
poles  in  the  first  quadrant  of  theV-plane  (cf.  Figure  14), 
and  the  arc  lengths  SLiT  and  S^T  are  shown  in  Figure  17. 

Thus,  the  a**1  longitudinal  or  transverse  Whispering  Gallery 

_ n 

wave  is  excited  at  the  critical  angle  &x  given  by 


(72) 


creeps  clockwise  ( SL^ ) or  counterclockwise  ( ) around 

the  cylinder  with  speed  otj 4 sine£T  and  radiates  off  at  the 
same  angle  after  n circumnavigations.  If  we  assume  that  the 
cylinder  is  slightly  lossy,  so  that  has  a small  imaginary 
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part,  then  as  a tends  toward  infinity,  SL>T  also  becomes 
infinite,  and  the  waves  which  creep  around  the  shadow  side 
of  the  cylinder  [including  the  multiple  circumnavigations 
( m*o)  ] are  exponentially  damped  out  and  never  reach  the 
observer.  Tne  waves  which  remain  are  those  which  traverse 
the  finite  arc  length  Su,t> 

-J  n»i  • 


(73) 


We  now  approximate  the  residue  sum  in  Eq.  (73)  by  an 
xntegral:  . 


(74) 


where  we  have  used  the  approximation2  4 for  ^ that  holds 
*or  larce  rv’. 


(75) 


The  approximation  of  the  sum  by  an  integral  is  justified 
because  goes  to  infinity  as  a tends  toward  infinity  and 

because  the  function  in  the  sum  oscillates  less  rapidly  as 
h increases,  so  that  the  contributions  to  the  sum  (integral) 
for  small  n tend  to  cancel  out.  Thus,  the  primary  contribution 
to  the  residue  sum  (integral)  comes  from  the  higher  order  poles. 
We  are,  therefore,  also  completely  justified  in  using  the  pole 
positions  of  Eqs.  (39) , since  the  method  used  to  calcu- 

late them  imposed  an  upper  limit  on  K which  increased  monot- 
onically  with  mode  number  (cf . Chapter  ID2.  and  Figure  5) . 
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We  point  out  that  the  method  of  approximating  the  residue  sum 
by  an  integral  is  similar  to  that  used  by  Tamir  and  Felsen25 
for  the  dielectric  slab  problem.  Rulf,22  on  the  other  hand, 
in  considering  the  fluid-fluid  (with  ) cylinder  problem, 

did  not  use  the  explicit  pole  positions  to  evaluate  the 
residue  sun,  but  converted  it  back  to  a contour  integral 
surrounding  the  poles.  We  treat  the  fluid-fluid  case  using 
our  method  in  Appendix  C. 

In  order  to  evaluate  the  integral  in  Eq.  (7  4)  , first  we 
change  variables 


and  then  we  evaluate  the  resultant  integral: 


(76b) 


where  we  can  assume  that  kL>r  has  a small  imaginary 
part  to  assure  convergence.26 

In  the  limit  of  infinite  radius,  the  expressions  for 
and  Pj.t  then  correspond  exactly  to  the  expressions 
for  the  longitudinal  and  transverse  lateral  waves  for  a flat 
elastic  half-space  [cf.  Eqs.  (B9) ] : 


P4>L  m 0-<> (p** -2°^) 


(77a) 
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P*>T 

with 

m 


i + kT  5 "0 

w 

(kr  Sr)* 

(77b) 


Conclusions 


We  have  established  the  connection  between  creeping  wave 
and  flat  surface  wave  theory  by  investigating  the  limit  of 
acoustic  scattering  from  a solid  elastic  cylinder,  imbedded  in 
a fluid,  whose  radius  tends  to  infinity. 

First,  we  calculated  the  asymptotic  behavior  of  the  com- 
plex circumferential  wave  numbers  by  substituting  the  appro- 
priate Debye-or  Airy-type  asymptotic  expansions  into  the  3x3 
secular  determinant  and  solving  it  using  iterative  techniques. 
The  creeping  wave  modes  fall  into  two  classes : those  with 

speeds  close  to  the  sound  speed  in  the  fluid  (Stoneley  and 
Franz  waves)  and  those  with  speeds  close  to  the  bulk  wave 
speeds  in  the  solid  (Rayleigh  and  Whispering  Gallery  waves) . 

It  was  found  that,  in  the  limit  of  infinite  cylinder  radius, 
the  wave  numbers  of  the  Rayleigh  and  Stoneley  modes  tend  toward 
those  of  the  Rayleigh  and  Stoneley  waves  on  a flat  elastic 
half-space,  while  the  Franz  mode  wave  numbers  tend  toward  the 
wave  number  of  sound  in  the  fluid.  The  longitudinal  and  trans- 
verse Whispering  Gallery  node  wave  numbers  tend  toward  the  long- 
itudinal and  transverse  wave  numbers  in  the  solid.  Graphical 
results  were  presented  for  an  aluminum  cylinder  in  water  (and 
in  one  case,  also  in  vacuum)  in  the  form  of  trajectories  in 
the  complex  wave  number  plane,  phase  velocities,  and  attenua- 
tions, all  as  functions  of  fluid  wave  number  times  cylinder 
radius.  The  results  show  good  agreement  with  existing  numerical 
results . 
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Then,  using  the  Watson-Sommerfeld  transformation,  we 
investigated  the  limiting  behavior  of  the  solution  to  the 
problem  of  the  scattering  of  a cylindrical  wave  from  a cylinder 
whose  radius  tends  to  infinity.  Using  the  analytic  expressions 
for  the  creeping  wave  numbers,  we  calculated  the  asymptotic 
behavior  of  the  residue  sums  corresponding  to  the  different 
classes  of  circumferential  waves.  It  was  found  that,  in  the 
limit  of  infinite  cylinder  radius,  the  Rayleigh  wave  for  the 
cylinder  goes  over  to  the  Rayleigh  wave  on  the  flat  elastic 
half-space,  while  the  Franz  and  Stoneley  waves  are  exponentially 
damped  out  (the  Stoneley  wave  contributes  to  the  field  in  the 
flat  case  only  at  glancing  incidence,  which  is  a special  case 
mathematically,  and  was  not  discussed) . In  the  limit,  the 
longitudinal  and  transverse  Whispering  Gallery  waves  combine 
to  form  the  longitudinal  and  transverse  lateral  waves,  respec- 
tively, for  the  flat  elastic  half-space. 

Thus,  the  transition  of  creeping  wave  to  surface  wave 
theory,  as  the  scattering  object  tends  toward  a flat  surface, 
has  been  established. 


Appendix  A 


Asymptotic  Expansions  of  Cylinder  Functions 


1 3 / 1 / 1 7 > 2 3 


In  the  following,  we  present  asymptotic  expansions  of 
cylinder  functions  where  both  the  values  of  the  (real)  argu- 
ment, -tj.,  and  of  the  modulus  of  the  (complex)  index,  V,  are 
large . 

1.  Debye  asymptotic  expansions 
Debye  expansions  are  appropriate  for  large  values  of 
= and  are  used  outside  the  circles  shown  in  Figure 

A1  whose  radii  are  determined  by  . The  Debye 


expansions 


Hf:\ V)  e*p  (L  ivcos-'v^-  ;«■/*} 

0 L ft/  , x »«/v 

Cf-y*-Vy 


^ $\f'-  e Xf>{- l CJf-yX) cos-'y/^  +L%j 
d W 


(Ala) 


x j1  + [_*  ai  (vi^) 


(Alb) 


are  valid  in  Regions  I and  II  of  Figure  Al,  with 

Jav'cj  C^-vV"(  < IT/a. 

T^e(cos",y/jN)£/n' 

Iwi  (cos *' v/^f')  4 O 0^ ''  <rv'- 

> O (.R.<t^'i<nviL)  . 
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(Ale) 
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These  regions  are  separated  by  the  curves 


(A2) 


and  by  the  portions  of  the  real  axis  as  indicated  in  the 
figure.  The  roots  of  and  of  lie  on  the  curves 

labeled  in  Figure  Al,  those  of  and  on 

the  curves  labeled  . 

In  Region  III,  one  has  instead 


U?(i)  ~ - i (-  <? 

..ft  1 r.  5_sv^r\  AA2 
x\_  1/ 


with 


WV'" 

(y*.  4 ft/x. 

^ (coslf'v/^  > O 
\X vw  (cosU~V^\  . 


(A3b) 


(A3c) 

In  Regions  IV  and  V,  the  appropriate  results  are  obtained 
by  using 

R.l^W'XVg) 

(A4) 

The  curves  h*t  are  the  anti-Stokes  lines  for  the 
asymptotic  expansions  of  Hy°^  , and  the  curves  h±x  those 
for  Hy(£  . in  the  vicinity  of  these  lines  (the  "transition 
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region") , the  sum  of  the  asymptotic  expansions  for  the  two 
adjacent  regions  has  to  be  used.  In  the  ensuing  "transition 
Debye  forms,"  the  two  exponentials  then  combine  to  form  trig- 
onometric functions,  while  outside  the  transition  regions,  one 
of  the  two  exponentials  is  dominant  and  the  other  subdominant. 

Debye  asymptotic  expansions  for  the  Bessel  functions  are, 
for  Region  III 


with 

\ 4 'fiy-a. 

Re(c > o 

\lw\(cos • (A5b) 

In  all  other  regions,  the  appropriate  expansions  are  found 
by  using 


together  with  the  previously  stated  results  for 

The  zeros  of  are  located  on  the  real  axis  in  Figure  Al, 


to  the  left  of 


1 • 


This  portion  of  the  real  axis  forms  the 


anti-Stokes  line  for  and  in  the  surrounding  transition 

region,  one  must  use 


X S{y\U) 
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with 

W=  (Jl-y^-yc os-'v/y'Zh^ 

where 

\om  i'S'/i. 

Re  (cos-'  v/y}  £ 4?  . 


(A7b) 


(A7c) 


2 . Airy-type  asymptotic  expansions 
For  the  case  that  , asymptotic  expansions 

of  the  cylinder  functions  which  are  expressed  in  terms  of 
the  Airy  function  are  appropriate.  They  are  valid  inside  the 
circles  of  Figure  Al.  We  have13'14 


+ fr|x' 

Hf - a^/y  c-  ^ei’FP)+ 


(A8a) 


(A8b) 


with 


and 


(A8c) 

(A8d) 

(A8e) 

(A8f) 

(A8g) 
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and  where  the  Airy  function  is  given  by 

C-  3.  ?3/z-  • 

^ 3 * > 

it  satisfies  Airy's  equation. 

For  Tvty , one  finds  the  compact  expression 


(A9a) 

(A9b) 


(A10) 


(All) 


Appendix  B 

The  Reflection  of  a Cylindrical  Wave  at  a Plane 


Fluid-Solid  Interface 


In  this  section,  we  calculate  the  total  acoustic  field 
due  to  a line  source  in  the  presence  of  a plane  fluid-solid 
interface,  where  both  source  and  receiver  are  situated  in 
the  fluid.  The  analysis  follows  that  of  Brekhovskikh2  and 


for  a point  source,  respectively.  The  field  can  be  resolved 
into  the  incident  wave,  the  geometrically  reflected  wave,  and 
surface  and  lateral  waves  which  propagate  along  the  boundary 
and  radiate  into  the  fluid. 

The  geometry  of  the  problem  is  shown  in  Figure  Bl.  An 
infinite  line  source  of  unit  strength  is  located  at  S.  The 
incident  pressure  wave  at  the  observation  point  P is,  there- 
fore®1 [with  a tine  factor  exp(-iajf)  suppressed] 


Using  the  Sommerfeld  representation1 2 for  the  Hankel  function, 
we  can  write  the  incident  field  as  a decomposition  into  plane 


Uberall,1  who  considered  the  fluid-fluid  and  fluid-solid  cases 


(Bl) 


waves  with  angle  of  incidence©": 

if t(*  so*  9-  +-  /i-iof 


(B2) 


where  the  contour C is  drawn  in  Figure  B2.  The  reflected 
wave  is  then  given  by1'2 


(B3a) 
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where  RCfi)  is  the  plane  wave  reflection  coefficient  for  a 
plane  fluid-solid  interface:1 


£=£-/iv  > 


(B3b) 


with  _ 

Z>+  = CO  s & [&*  - S iV'G)  + s ' 

and 

“ 6<sr" bSr/h.  • 


(B3c) 


(B3d) 


Equation  (B3a)  for  the  reflected  wave  can  be  rewritten 
as 

J=Wfi  = dfr  R<»  eifcf 

c (B4) 

where  is  the  distance  from  the  image  source  S1  to  the 
observation  point  P. 

The  integral  in  Eq.  (B4)  is  evaluated  using  the  method 
of  steepest  descent1'2,  where  we  assume  /?|0/>>1  . The  path 
of  steepest  descent  C5  and  the  saddle  point  (corresponding 
to  the  angle  of  incidence  or  reflection  of  the  geometrically 
reflected  wave)  are  shown  in  Figure  B2.  The  result  for  the 
geometrically  reflected  wave  is  then 


W 


(B5) 


as  expected  from  geometrical  acoustics. 

Additional  contributions  to  the  reflected  field  arise 
from  any  singularities  of  R{&)  which  are  crossed  as  contour 


55 


C is  deformed  into  contour  C3.  Branch  point  singularities, 
provided  by  the  radicals  r are  found  at 


- ± s i ^ 0-T-  1 Si*-'  «T 


(B6) 


while  poles  [corresponding  to  solutions  of  D+=0,  which  is 
completely  equivalent  to  the  generalized  Rayleigh  Eq.  (8)] 
of  physical  interest  occur  at 

&*r  s,vr‘  3 &s  - s'V(ks(  k)  (B?) 

where  and  It*  are  the  Rayleigh  and  Stoneley  wave  numbers 
[e.g.,  Eqs.  (11)  and  (16)  for  the  water-aluminum  interface). 

These  branch  points  (and  the  corresponding  branch  cuts)  and 
poles  are  shown  in  Figure  B2 . 

The  angles  Re  9^,  Rt©r,  ©t,  and  are  the  critical  angles 
of  incidence  for  excitation  of  the  Rayleigh,  Stoneley,  and 
longitudinal  and  transverse  lateral  waves,  respectively  (for 
the  water-aluminum  interface:  ©l_=13°26’,  ©j.=29°20  ' , R*  £^=31°36  ' , 

Re  6^=90°)  . When  dl,  and  frT , as  in  Figure  B2,  then  the 

Rayleigh  and  lateral  waves  contribute  to  the  field  at  the 
observation  point  P.  The  Rayleigh  wave  contribution  is  found 
from  the  residue  at  the  Rayleigh  pole,  while  the  longitudinal 
and  transverse  lateral  wave  contributions  arise  from  integrals 
around  the  branch  cuts  originating  at  and  <9^-,  respectively. 

The  resulting  expressions  are: 

Rayleigh  Wave: 


(B8a) 
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where 

\-o~  *o|cos%_  3 L-  */cos^  i>+(&)-  (B8b) 


Longitudinal  Lateral  Wave: 

e-#  &i[hO*i-0+iLLj 

Pl  7*1*  ^ 


(B9a) 


Transverse  Lateral  Wave: 

? = 3 |s';^  ei&ft.<W-feTLT] 

rT  **  <Ur)3/*  ““ 

(B9b) 

where 

L0-  ^o/cosB*ST  ^ L-^/cos8i>T  > m = /Vo  * (B9c) 

The  geometrical  meaning  of  these  results  is  clear  from 
inspection  of  the  phases  and  Figures  B3  and  B4.  In  each 
case,  the  wave  is  excited  at  its  critical  angle,  propagates 
along  the  interface,  and  reaches  the  observation  point  P by 
radiating  into  the  fluid  at  the  same  angle. 

The  Stoneley  wave  is  not  excited  unless  6^**®/*  (glancing 
incidence) , in  which  case  the  method  of  steepest  descent  must 
be  modified  to  take  into  account  the  effect  of  a pole  near 
the  saddle  point.2  This  case  will  not  be  discussed  here. 


Appendix  C 

Creeping  Waves  and  Lateral  leaves 
for  the  Fluid  Cylinder 


In  the  case  of  scattering  by  a fluid  cylinder  (CT-»0)  we 
have  the  residue  sun  "p*.  of  Eqs.  (53)  or  (61)  with 

CL^- 

p-  (CD 

The  Whispering  Gallery  pole  positions  are  given  by14 

VU)L^=  + fr(f)  (c2) 

where,  as  in  the  case  of  the  solid,  the  material  properties 
enter  only  in  higher  order  terms.  Using  the  appropriate 
Debye  expansions  for  the  Hankel  functions  of  X and  Airy-type 
expansions  for  the  Bessel  functions  of  (cf.  Appendix  A) 
along  with  Eqs.  (65)  and  the  fact  that 

?y  \ — j l°a-  /i..«va\l/x 

v 


(C3) 


we  obtain  the  result 

l-l 


ssj 


- 


^ (C4) 

which  is  just  the  limit  of  Eq.  (67)  for  the  solid  as  o(T 
goes  to  infinity.  Then,  using  Eqs.  (62),  (46b),  (70),  and 
(C4) , and  approximating  the  residue  sum  by  an  integral  as 


before,  we  find  that 


i •+■ fei.sC] 


?1,L  'JW  mQ -«i) 


(C5) 
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which  is  just  the  limit  of  Eq.  (77a)  as  goes  to  infinity 
and  corresponds  exactly  to  the  expression  for  the  lateral 
wave  for  a flat  fluid  half-space2 ' 2 2 . 


List  of  Illustrations 


Complex  V -plane  with  various  regions  in  which  dif- 
ferent asymptotic  expansions  of  cylinder  functions 
are  employed  in  our  analytic  calculation. 

Rayleigh  zero  trajectory  in  the  complex  plane 
at  varying  values  of  the  parameter  ka,  for  the  case 
of  an  aluminum  cylinder  in  water. 

Crosses:  present  work;  solid  circles:  numerical 

results  of  Ugincius9;  solid  square:  flat  Rayleigh 

limit.  Also  shown  is  the  location  of  hrfk.  • whose 
(real)  numerical  value  is  0.491. 

Stoneley  zero  trajectory  in  the  complex  plane 
at  varying  values  of  the  parameter  ka,  for  the 
case  of  an  aluminum  cylinder  in  water. 

Franz  zero  trajectories  in  the  complex  plane  kprxjfe.  at 
varying  values  of  the  parameter  ka,  for  the  cases  of 
soft  (■) , rigid  (•) , and  aluminum  cylinders  (x)  in 
water:  (a)  first  Franz  zero,  v\=l,  (b)  second  Franz 

zero,  A=2 . 

The  value  of  ka,  as  a function  of  mode  number  n,  at 
which  the  absolute  value  of  the  relative  error  in 
the  approximation  to  the  tangent  function,  used  to 
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Fig.  6 


Fig.  7 


Fig.  8 


Fig.  9 


calculate  the  longitudinal  Whispering  Gallery  zeros, 
begins  to  exceed  25  percent. 

Rayleigh  zero  trajectory  in  the  complex  plane 
at  varying  values  of  the  parameter  ka,  for  the  case 
of  an  aluminum  cylinder  in  water.  Crosses:  analytic 

results;  square:  flat  Rayleigh  limit  (ka-»«>)  . Circles 

numerical  results. 

(a)  Normalized  phase  velocity  ^'/c.  of  the  Rayleigh 

wave  for  an  aluminum  cylinder  in  water,  plotted  vs. 
ka.  Circles:  Numerical  results,  curve:  analytic 

results. 

(b)  Normalized  phase  velocity  of  the  Rayleigh 

wave,  plotted  vs.  k*a  . Circles:  numerical  results 

for  aluminum  cylinder  in  air.  Curve:  analytic 

results,  and  crosses:  Viktorov' s 8 ' 1 8 results,  both 

for  aluminum  cylinder  in  a vacuum. 

Dispersion  curve  for  the  Stoneley  wave  on  an  aluminum 
cylinder  in  water.  Solid  circles:  numerical  results. 

Crosses:  analytic  results. 

Stoneley  and  first  five  Franz  zeros  for  an  aluminum 
cylinder  in  water.  (a)  Dispersion  curves,  (b)  normal- 
ized attenuation.  For  the  Franz  zeros,  analytic 
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results  are  given  by  solid  curves,  numerical  results 
by  open  circles.  For  the  Stoneley  zero,  analytic 
results  are  given  by  crosses,  numerical  results  by 
solid  circles. 

Fig.  10  (a)  Dispersion  curves  for  the  first  five  transverse 

Whispering  Gallery  zeros  and  the  Rayleigh  zero  for  an 
aluminum  cylinder  in  water.  The  numerical  results  are 
given  by  solid  curves;  the  analytic  results  are  shown 
as  long  dashes  where  they  are  most  accurate  and  short 
dashes  where  they  are  less  reliable. 

(b)  Normalized  attenuations  for  the  Whispering  Gallery 
zeros.  Solid  curves:  numerical  results  for  first  four 

zeros.  Short  dashes:  extrapolated  numerical  results. 

Long  dashes:  analytic  results  for  first  two  zeros. 

Fig.  11  The  first  seven  longitudinal  Whispering  Gallery  zeros 
for  an  aluminum  cylinder  in  water.  (a)  Dispersion 
curves,  (b)  normalized  attenuation.  The  analytic 
results  are  shown  as  solid  curves  where  they  are  most 
accurate  and  dashes  where  they  are  less  reliable. 

Fig.  12  Geometry  of  the  cylinder  scattering  problem;  line 
source  at  S,  observer  at  P. 


Fig.  13  The  complex  V -plane  showing  contours  for  the  Watson- 

Sommerfeld  transformation  used  in  the  elastic  cylinder 
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Fig. 


Fig. 


Fig. 


Fig. 


Fig. 


scattering  problem  [Figure  taken  from  Doolittle, 
et  al.1*  with  permission  of  the  American  Institute 
of  Physics] . 

14  The  complex  V-plane  showing  contours  for  separating 
out  the  geometrical  wave.  Also  shown  schematically 
are  the  positions  of  the  saddle  point  and  the 
Rayleigh  ( fe.^0.)  , Stoneley  ( kj»d-)  , Franz  (vF)  , and 
longitudinal  ( ) and  transverse  ( V,jT)  Whispering 
Gallery  poles  for  a large,  fixed  value  of  ka. 

15  Franz  or  Stoneley  wave  getting  on  the  cylinder  tan- 
gentially, creeping  clockwise  (Fig.  15a)  or  counter- 
clockwise (Figure  15b)  around  it,  and  getting  off 
the  cylinder  tangentially. 

16  Rayleigh  wave  being  excited  at  the  critical  angle 
creeping  clockwise  or  counterclockwise  around  the 
cylinder,  and  radiating  off  at  the  same  angle. 

17  Longitudinal  or  transverse  Whispering  Gallery  wave 
being  excited  at  the  critical  angle  <9^r  , creeping 
clockwise  or  counterclockwise  around  the  cylinder, 
and  radiating  off  at  the  same  angle. 

Al  Complex  V -plane  showing  regions  of  validity  for  dif- 
ferent asymptotic  expansions  of  the  cylinder  functions, 
used  in  the  analytic  calculation. 
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Fig. 


Fig. 


Fig. 


Fig. 


B1  Geometry  of  an  observer  at  P receiving  a cylindrical 
wave  from  a line  source  at  S and  a reflected  wave 
from  the  image  source  at  S'. 

B2  Integration  path  C for  the  incident  and  reflected 
waves  in  the  complex  £”-plane;  saddle  point  6i  with 
path  of  steepest  descent  Cg;  Rayleigh  pole  , Stoneley 
pole  , and  branch  points  ©<.,  0-r  with  corresponding 
branch  cuts  (dashed  lines) . 

B3  Cylindrical  wave  from  source  S exciting  Rayleigh  wave 
at  point  A,  which  propagates  along  interface  and 
reaches  observation  point  P by  radiating  into  fluid 
at  point  B. 

B4  Cylindrical  wave  from  source  S exciting  longitudinal 

or  transverse  lateral  wave  at  point  A,  which  propagates 
along  interface  and  reaches  observation  point  P by 
radiating  into  fluid  at  point  B. 


Fig. 


0.030 


0.025  - 

k 

0.020  - 

0.015  _ 
0.43 


ko  = l5 

ka  = 40 

x 


25 


I 


1 


0.44  0.45  0.46  0.47 


50 

X 


40 


1 


60 

x 


70 

*80 
x 90 


50 


60 


xxIOO 
70  80 


1 


200 
x 300 


RAYLEIGH 

LIMIT 


1000 


1 


500 


1 


0.48 

kRV 


0.50  0.51  0.52  0.53 


cry 

Ln 


H- 

*0 


k 


ro 


\ 


0_x  2000 


JP 


x 3000 

4|-x  4000 

X 5000 

x 6000 
2 -X  7000 

9000  SS  8000 
9000  x 10,000 

^STONELEY  LIMIT 

1.00  1.02  1.04 


90 

XXX  X 
100  80  70 


_J I 

1.06  1.08 


i 

1.10 

kS'r/k 


X 

60 


_L_ 

1.12 


x 

50 


L_ 

1.14 


x 

ka«40 


_J I I 

1.16  1.18  1.20 


<T\ 


U) 


Fig. 


67 


Fig.  4a 
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Fig.  4b 

0.20 

1— 

0.18 

F2 

■ ka 

= 60 

0.16 

— 

80 

60#" 

v 

■ 70 

014 

A 

70 

• 

v 

0.12 

80 
90  * 

ill  o 

& |jg 

0.10 

0 

,0°*  X 

0.08 

— 

200 

a 

0.06 

0.04 

300 

* 

400 

500* 

0.02 

1000 

- * 

10,000 

* 1 1 

_J L_ 

1 

1 1 

1.00  1.02  1.04 

1.06  1.08 

1.10 

1.12  1.14 

v_  Fr 
ka 

1200r 


1000 


800 


cr\ 

vo 


k0MAX  600 


400 


200 


0 


5 6 7 

MODE  NUMBER  n 


8 


10 


11 


Ln 


Fig. 


o 


k 


CT\ 


Fig. 


— 


— 


— 





_ 


Fig.  7a 


72 


Fig.  7b 
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Fig.  17 
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Fig.  B4 
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